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FORMATION AND PROPAGATION OF A SHOCK WAVE 
FROM A CONCAVE BOUNDARY 
by Dudley G. McConnell 
National Aeronautics and Space Administration 

SUMMARY 

This report shows when and where a shock wave forms in a two-dimensional flow 
field and presents a method for predicting the shape and strength of the shock wave. 
Knowledge of the point of formation, the strength, and the shape of shock waves are im- 
portant in the design of inlets and the design of blading for turbomachines. 

The present method consists of first finding the envelope of characteristics forming 
in the flow field from the bounding streamline, as suggested by Tollmien. In a simple 
wave field, where one family of characteristics consists of straight lines, the envelope 
can be given in closed form; in a complex wave field, where both families are generally 
curved, the envelope must result from a solution of the whole flow field. Once the 
envelope is found, criteria and methods are presented for locating the shock within the 
envelope. In simple waves, the shock path results from a point by point integration; in 
complex waves, the principle of minimum entropy production may be used to pick out 
the most likely shock location. Specific examples of the simple wave case are pre- 
sented. 


INTRODUCTION 

The object of this report is to solve the following problem: For given upstream 
conditions and a given boundary profile, determine whether a shock wave will occur and, 
if so, find the location, strength, and shape of the shock wave. 

The motivation for this study comes from the current interest in supersonic inlets 
and supersonic operation of turbines and turbocompressors. Focused characteristic 
inlets have been successfully designed for some time; however, at high supersonic 
Mach number hidden shocks appear. These unexpected shock waves, which decrease 


the efficiency of the inlet, have not yet been analyzed. Furthermore, there is interest 
in finding the path of shock waves that may occur during off-design operation. There 
is a similar need to study the formation of shock waves within turbomachinery blade 
rows and the propagation of shocks between stages. Impingement of a shock wave on 
a blade could radically alter blade loading and flow incidence. This off-design loading 
may enhance blade flutter and lead to early blade failure. Off-design flow incident 
would impair performance. 

M. B. Abbott (ref. 1) shows that shock-wave formation always results from the 
convergence of characteristics. As a compression wave processes the gas, the speed 
of sound increases behind the wave. In a train of compression waves, for instance, those 
coming from a concave boundary or an accelerating pistion, trailing waves tend to over- 
take the leading wave. If these trailing waves coalesce into the lead wave, the wave 
front steepens, and its speed increases. As coalesecence continues, the process through 
the wave finally becomes discontinuous; that is, a shock wave forms. 

In steady, plane, supersonic flow, characteristics are propagation paths for Mach 
waves. So the coalescence of trailing waves shows up graphically as the formation of 
an envelope of characteristics (refs. 2 and 3). This process is shown in figure 1. 

Usually the envelope will be cusped in its initial point. The initial point of the envelope, 
the cusp, is also the initial point of the shock wave. Courant and Friedrichs (ref. 2) 
Oswatitsch and Kuerti (ref. 3) and others discuss the general nature of the envelope but 
do not give details on its calculation. Tollmien (ref. 4) gives more detail but then con- 
centrates on a series extension of the analytical solution downstream of the envelope. 

He does not calculate the shock shape within the envelope. 

The present treatment considers a characteristic envelope and the resulting shock 
wave that form in the external flow over a concave profile. The envelope forms in a 



Figure 1. - Formation of envelope in simple wave flow. 
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simple wave region, so its shape and its starting point can be given in closed form. A 
point by point integration fixes the shock path once the initial point of the envelope is 
known . 

SYMBOLS 

A mean slope along shock wave, defined after eq. (7a) 

B slope of characteristic, defined after eq. (7a) 

EpEg parametric form of locus of characteristic envelope 
F function giving location of bounding streamline 

G parametric form of family of characteristics 

K j/(y + WW ~ 1) 

M Mach number 

m slope 

q magnitude of flow velocity 

R channel radii of curvature^ 

s entropy 

x, y spatial coordinates, physical plane 1 

a Mach angle 

y ratio of specific heats 

rj characteristic net variables 

Ti pressure ratio across shock wave 

i? flow deflection angle 

co Prandtl-Meyer function 

Subscripts: 

A upstream border of characteristic envelope 

B return flow region of characteristic envelope 

C downstream zone of characteristic envelope 

e envelope 

^In analysis of channel flow, spatial distances are measured in entrance half 
heights . 


3 



p boundary profile 

00 condition of undisturbed flow 


ANALYSIS 

Envelope Conditions and Shape 


If, as in figure 1, a compression region borders a region of constant state, at least 
one family of characteristics will consist of straight lines. This family of lines can be 
written as 


G(x,y;x p ) = 


- F <V - 


(x 


i 


tan(t? 


+ cr 


= 0 


(1) 


where F(x p ) is the contour of the bounding profile or bounding streamline. For given 
upstream conditions, the angle or of each Mach line to the approaching flow direction 
is a function of profile position x through 

ir 



If an envelope exists, then neighboring characteristics will intersect along the envelope. 
And a characteristic will intersect the envelope at a point of tangency. According to 
Courant (ref. 5), these conditions are both necessary and sufficient, subject to the con- 
straint that both the envelope and the members of the family in fact have well defined 
tangents. Thus, since a well defined characteristic field is given, the conditions for 
the existence of an envelope are in fact conditions on the shape of the bounding stream- 
line. 


Assume that an envelope exists and is given by 

x e = E 1<V 

= E 2<V j 


( 2 ) 


Some characteristic intersects the envelope in the point (x°,y°) such that (x°,y°) satis- 
fies the equation 


-( : »p) = 0 
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So we may write 


3G/ dE l\ + aG/ d 52 \ + /0G 1= 0 


ax \dx p y ay yix p/ 


1 3x„ 


The slope of the characteristic through the point of intersection is 


dy = . Qo/3*) 

dx (3G/ay) 


The slope of the envelope is 

/dyX (dE 2 /dx p ) 
\ dx e/ ( dE i/ dx p ) 


If this is in fact an envelope; these slopes must be equal. Consequently, 


9G 

0X 




= 0 


(3) 


(4) 


(4a) 


The result is that 


0G 

3x„ 


1= 0 


(5) 


which along with G(x,y;x ) = 0 presents two equations to relate x and y to x along 

r r 

the envelope. 

For the characteristic family from equation (1), the envelope is generally given by 


<*e - x p ) = 


tan(^ p + a p ) - tan ,? p 




sec (^ p + a p ) 


da 

1 +( —2 
ids 


/di? 
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(y e - y p ) = 


|tan(i? p + a p ) - tan ^ p Jtan(^ p + a p ) 
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Notice that the curvature (d^p/dXp) must be positive (i. e. , a concave profile) for an en- 
velope to exist in a supersonic compression flow. If the curvature is infiinte, as in a 
corner, the envelope begins on the profile itself. 


Shock Wave Region 

Tollmien (ref. 4) shows that the shock wave and the envelope have a common starting 
point: the point for which, in this case, 



= 0 


Given this initial point, it is possible to integrate the equation 



(6a) 


(7) 


point by point along the shock wave - j being a known function of the upstream Mach 

s 

number and the flow angle 

An alternate difference technique could use a simultaneous solution of the equations 

ys,m-l-y s ,n _ tan Wl atan9 s,n'l 
x s , n+1 " X s, n 

^s,n+l ^p,n+l 
x s,n+l x p,n+l 

which amounts to about the same thing. 

Figure 2 shows the envelope and shock wave which develop in a flow that approaches 

/ -axp\ 

a profile of the form F(x_) = 1 - e * up to its point of inflexion. This curve has an 
2 P \ / 

inflexion at x^ = l/2a; beyond this poim the profile is continued as a straight line at the 
angle 


tan(,9 i+O! . ) 
' p, n+1 p,n+r 


(7a) 
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(a) Approaching Mach flow, 2.0. (c) Approaching Mach flow, 6.0. 

Figure 2. - Comparison of envelope and shock wave paths. 


For the conditions given in figure 2, the shock and the envelope agree closely up to flow 
deflections of about 5° to 8°. This is equivalent to the region where the Prandtl-Meyer 
epicycloid approximates the shock polar. As a measure of the accuracy of the shock 
path, the far field shock angle, using the present method, was within 1° of the tabulated 
shock angle for two-dimensional flow (ref. 6) for the range of Mach numbers considered 
The main contribution of this section has been the integration for the shock path 
using the initial point of the envelope as the starting point of the integration. 


Complex Wave Region 

Simple wave shock formation occurs most often in external flows. Internal flows 
usually lead to complex wave interactions since all boundaries affect the flow. However, 
in a complex wave field, neither family of characteristics is a set of straight lines. 

Thus, the existence of an envelope can only be known from a complete solution of the 
flow problem, that is, from a solution of the fundamental characteristic equations 
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( 8 ) 



Along with appropriate initial or boundary data, these equations furnish the functions 
x(|, rf) and y(£, 77) for mapping the flow variable q onto the physical plane 


Characteristic net Physical plane"'* 


q(£, vf 

► 


q(x, y) 

4x,y) 


> 


j 


(9) 


The techniques for finding the mapping functions are given in references 1 to 4. There- 
fore, assume that such a mapping exists. And further, following Tollmien (ref. 4), 
assume that the envelope forms because of the convergence of left running characteris- 
tics (77 = constant). (See fig. 3.) Along such a characteristic, £ is the current variable. 
Each such characteristic will touch the envelope in one point £ e = ^( 77 ). Then the 
mapping of the envelope in the physical plane is 

x e = x[s e (7?),'/7] 
y e = y[y ^ 


Where each characteristic touches the envelope, the slope of the characteristic will equal 
the slope of the envelope. Thus, the slope along the characteristic is 

/dy\ = (dy/dj) 

\ dx / ??= const 

This must equal the slope of the envelope at that point 

dy„ (3y/9^)(dA/d77) + (?y/3 77) 

— - = ? ( 10 ) 

dx g (3x/3|)(d| e /d7j) + (3x/3 77 ) 
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Cross multiply to obtain 



= 0 


This equation may be written as 



(10a) 


Exclude the exceptional case (9x/3£) = 0, which simply implies vertical characteristics. 
Then, the general condition for the envelope is 



( 11 ) 


This condition shows the steepening process at the wave front. It states that two Mach 
waves coincide at one (x,y) position. 

For cases of interest, the envelope will form within the flow field. Thus, there 
will be some extreme value beyond which Z e (v) will not go (fig- 3)- Whether this ex- 
treme value is a minimum or maximum (as in ref. 4) depends on the curvature of the 
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boundary generating the wave train. In any event, at the extreme value, following ref- 
erence 4, a sufficient condition is 



( 12 ) 


Let the boundary of the envelope be given by 


— = nuv),v) = 0 

d V 


Then 


d V \3 J VvJ 

So at the extreme value, the initial point 


■\ 


dF 

d77 


_ 9 x _ q 

dv 2 


3 x 

a , 3 


* 0 


(13) 


(This result excludes the unlikely case in which (9x/3 rj) becomes infinite. ) Now the loca- 
tion of the envelope is known in both simple and complex wave regions. But within the 
envelope region the path of the shock wave is still unknown. 

The final determination of the shock wave rests on the minimum entropy production 
theorem (ref. 5). The reasoning is as follows (see fig. 4). Within the envelope, the 
mapping x = x(|, rj), y = y(£, rj) is triple valued. The envelope region B would be a field 
of return flow. At any value of (x, £) there are three states available to the gas: 77 ^, 

77 g, 77 ^,. According to the second law, the shock process will take the gas to the avail- 
able state with the locally maximum entropy. Hence, from branch A to branch C. 
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Figure 4. - Typical plot used to determine shock wave path. 


Branch B is not allowed since it would imply a returning flow. But it is still necessary 
to choose one out of a continum of points such that x^ and x^, coincide. The minimum 
entropy production theorem will pick out the most likely value of x^ — x^, in which the 
shock will occur. 

Since the Mach angle a increases through a shock wave, (Stj/Dx^ > (drj/dx)^. 
Consequently, the jump in 77 gets larger at larger downstream positions x. But the 
jump in 77 is proportional to the flow deflection if £ is constant. 




C 



Since shock strength increases with flow deflection, minimum entropy production im- 
plies that the shock will occur at the earliest possible location. Hence, the shock wave 
follows the envelope. 

The entropy jump across the wave is found from 


AS 

R 



- In 77 
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where 



and 


«(Mj = — 

A 2 

As an example, consider the steady plane flow contained between two circular arc 
walls (fig. 5). The inlet conditions are ^ = 0 and = 1. 503. The characteristic 



Figure 5. - Physical layout of channel showing shock wave location. Note: arc 
BAC is the beginning of the complex wave field shown'in figure 6. Upstream 
flow conditions: Mach 1. 503; flow deflection angle, a 


map (a point plot) is also shown (fig. 6) as it was obtained from a machine solution and 
electronic plotter. The point plot clearly shows the convergence of the left running 
waves. This convergence occurs in spite of a reflection of the compression waves off 
the upper expansion surface. Indeed, the envelope forms in a generally expanding flow; 
the inlet Mach number is 1. 503, but the envelope forms at = 1.92 and = 25°. 
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Figure 6. - Computer point-plot showing characteristics 
converging. Upper arc, 12; lower arc, 16; Mach L5; 
flow deflection angle, 0, 


Figure 7 is a plot of the physical coordinate x as a function of the characteristic 
coordinates £ and 77 . The initial point of the envelope is the point where 



= 0 



- 0 


and 



which figure 7 verifies. For the local flow conditions of this example (M^ = 1.92) there 
is a very small entropy change in the vicinity of the initial point, (AS/R) - 10 “ 4. Never- 
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AXIAL COORDINATE, 





theless, the pressure change across the jump is not negligible (77 = 1. 15). Figure 8 shows 
the locus of the shock wave in detail. 

SUMMARY AND CONCLUSIONS 

This report presents methods for predicting shock formation and shock-wave tra- 
jectory in simple and complex wave fields. In simple wave regions, say for external 
flows, the method consists of point by point integration, given the shock slope by either 
method, starting out from the initial point of the characteristic envelope. In complex 
wave fields, minimum entropy production principle indicates that the shock wave lies 
along the upstream boundary of the characteristic envelope. 

At the present time these methods are being extended to technical applications. In- 
tegral boundary-layer solutions are being used (refs. 7 and 8), along with the simple 
wave method, to predict the shock trajectory that causes a laminar boundary layer to 
separate as it approaches a concave bend. The method for complex wave flows is being 
extended to include the effect of envelope interaction with a boundary and to include the 
effect of channel exit on the shock wave path. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 31, 1971, 

129-01. 
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